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We present an analytical model capable of describing two-photon ionization of atoms with attosec¬ 
ond pulses in the presence of intermediate and final isolated autoionizing states. The model is based 
on the finite-pulse formulation of second-order time-dependent perturbation theory. It approxi¬ 
mates the intermediate and final states with Fano’s theory for resonant continua, and it depends on 
a small set of atomic parameters that can either be obtained from separate ab initio calculations, or 
be extracted from few selected experiments. We use the model to compute the two-photon resonant 
photoelectron spectrum of helium below the N=2 threshold for the RABITT (Reconstruction of 
Attosecond Beating by Interference of Two-photon Transitions) pump-probe scheme, in which an 
XUV attosecond pulse train is used in association to a weak IR probe, obtaining results in quan¬ 
titative agreement with those from accurate ab initio simulations. In particular, we show that: i) 

Use of finite pulses results in a homogeneous red shift of the RABITT beating frequency, as well 
as a resonant modulation of the beating frequency in proximity of intermediate autoionizing states; 
ii) The phase of resonant two-photon amplitudes generally experiences a continuous excursion as a 
function of the intermediate detuning, with either zero or 27r overall variation. 

BAGS numbers: 32.80.Qk,32.80.Fb,32.80.Rm,32.80.Zb 


I. INTRODUCTION 

In the photoionization of poly-electronic systems, ab¬ 
sorption of an energetic photon is often associated to the 
formation of autoionizing states with a lifetime of few 
femtoseconds. Until recently, the role of such states in 
photoemission could only be studied in stationary con¬ 
ditions, typically using the long pulses (tens of picosec¬ 
onds) generated at synchrotron facilities [1]. Table-top 
attosecond sources [2-5], which deliver extreme ultravi¬ 
olet (XUV) pulses capable of coherently exciting several 
states in the continuum across wide energy ranges [6-8], 
associated to control pulses within pump-probe schemes, 
have opened the possibility of studying the dynamics of 
metastable wavepackets at its natural time scale. For 
example, it has been possible to follow in time the grad¬ 
ual depletion of individual autoionizing states [9] , as well 
as the progressive buildup of their population across the 
pump sequence [10], and to reconstruct the rapid evo¬ 
lution of an autoionizing wavepacket [11] from beatings 
between its individual metastable components. 

Among attosecond interferometric spectroscopies, a 
prominent role is occupied by the so-called Reconstruc¬ 
tion of Attosecond Beating by Interference of Two- 
photon Transitions technique (RABITT) [12, 13], which 
makes use of weak pump and probe pulses and is thus 
amenable to a perturbative treatment. In RABITT spec¬ 
troscopy a target atom or molecule is ionized by a train 
of attosecond pulses (APT), acting as a pump, in associa¬ 
tion with a weak long IR probe pulse, with a controllable 
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time delay r between APT and probe. The spectrum 
of the APT, which is generated from the interaction of 
an intense IR pulse with an active medium [12, 14], is 
formed by odd harmonics u) 2 n+i of the fundamental IR 
frequency, uj 2 n+i = (2n -I- 1 ) wir , while the IR probe is a 
weak replica of the IR pulse used to generate the train. 
When the APT pump and the IR probe overlap, there¬ 
fore, the target can either absorb one XUV photon from 
harmonic 2n — 1 plus one IR photon, or absorb one XUV 
photon from harmonic 2n -|- 1 and emit, in a stimulated 
way, one IR photon. These two processes interfere, giv¬ 
ing rise to a sideband photoelectron signal SB 2 n which, in 
stationary conditions, oscillates as a function of the time 
delay as Isb 2 „ = 4 cos(2a;iRT -|- A</)h -I- Atpat) [15], where 
A<()h is the phase difference between two consecutive har¬ 
monics in the APT spectrum, while the so-called 

atomic phase, is the relative argument of the two-photon 
transition matrix elements for the IR absorption and the 
IR emission quantum paths. 

If Atpat is a known slowly varying function of photo¬ 
electron energy, from the beating of the RABITT side¬ 
bands one can recover the relative phase between the 
harmonics in the train. Use of the RABITT technique 
with this approach has been instrumental to demonstrate 
that the harmonics from High-Harmonic Generation [16] 
(HHG) came in the form of trains of pulses [12], to under¬ 
stand the generation of attosecond light bursts [17] and 
to develop phase-compensation techniques that minimise 
the duration of individual pulses within the train [18, 19]. 
Gonversely, if the harmonic phases are known, from the 
sideband beatings it is possible to reconstruct the atomic 
phases [20]. This latter approach permits one to mea¬ 
sure both phase and amplitude of the dipole transi¬ 
tion matrix element from the ground to the intermedi- 
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ate continuum states and, in turn, to reconstruct the 
dynamics of the photoemission event. This procedure 
has been followed, for example, to determine the rela¬ 
tive delay between photoemission from the 3s and the 3p 
sub-shells of argon [21, 22], the phase difference between 
photoemission from the outermost s shell in different rare 
gases [23, 24], the energy-sharing resolved double ioniza¬ 
tion of Xenon [25] , and the nuclear dynamics in H 2 [26] . 
In all these examples, the intermediate continuum states 
do not feature any distinct structure. The RABITT tech¬ 
nique, therefore, could be used to infer properties of the 
continuum that vary smoothly across the energy gap 2a;iR 
separating two consecutive sidebands in the photoelec¬ 
tron spectrum. 

If, on the other hand, one of the two-photon paths in 
the RABITT scheme is nearly in resonance with a nar¬ 
row intermediate bound or metastable state, both the 
amplitude and phase of the corresponding sidebands ex¬ 
hibit a strong modulation as a function of the detun¬ 
ing of the energy of the harmonic closest to the reso¬ 
nance [27, 28]. When the contribution of the continuum 
to the resonant path is negligible, the sideband phase- 
shift undergoes a jump of tt as the frequency of the har¬ 
monics, increased gradually, traverses the energy of the 
intermediate resonant state. This phenomenon was ob¬ 
served experimentally with helium, using the ls3p Ryd¬ 
berg state as intermediate resonance [29], as well as in 
the N 2 molecule [30, 31], where the intermediate reso¬ 
nance was an autoionizing vibronic state. Even in the 
latter case, the contribution of the intermediate contin¬ 
uum turns out to be negligible, despite the fact that the 
resonance in N 2 does interact with the ionization channel 
to which it eventually decays. 

In general, however, both the continuum and the 
localized part of an intermediate resonant state may 
contribute to the two-photon transition. In a recent 
work [10], we showed that in such case, instead of un¬ 
dergoing a distinct jump of tt, the sideband phaseshift 
exhibits a finite excursion. For short pulse durations, 
furthermore, not only is the phase of the sideband os¬ 
cillations affected by an intermediate resonant state; the 
frequency of the sideband beating is modified as well. 
For very short pulses, in fact, even in absence of an in¬ 
termediate resonant state, the beating frequency is red 
shifted with respect to the 2a;iR nominal value. In the 
present work we provide a full derivation of an analyti¬ 
cal model, first presented in [10], which explains all these 
phenomena. The model provides resonant second-order 
ionization amplitudes, in which both the intermediate 
and the final continuum states of the two-photon transi¬ 
tion matrix element can feature isolated resonances, each 
described with Fano’s formalism, and the external field 
is formed by an arbitrary number of Gaussian pulses. 
In [10], by comparing the predictions of such model with 
accurate ab initio simulations [10, 32-34], we showed that 
this approach permits us to reproduce quantitatively the 
ionization of the helium atom to the region of the N = 2 
autoionizing states with a sequence of attosecond pulses 


in association with a weak IR probe pulse. 

The paper is organised as follows. In Sec. II we give 
an overview of the main formulas for two-photon tran¬ 
sition amplitudes with finite pulses, we introduce and 
justify the on-shell approximation for the calculation of 
continuum-continuum transition matrix elements, and 
comment several aspects of non-resonant RABITT with 
finite pulses, including an explanation of the red shift of 
the sideband beating. In Sec. Ill we derive the finite- 
pulse two-photon resonant model in the single-channel 
case as well as some straightforward generalisations. In 
section IV we apply the model to the resonant photoion¬ 
ization of helium with the RABITT technique and com¬ 
pare the analytical predictions with the numerical results 
we obtain with accurate ab-initio numerical simulations. 
In Section V we draw our conclusions. 


II. TWO-PHOTON TRANSITIONS 

In this section we derive the lowest-order perturbative 
expression for two-photon transition amplitudes with fi¬ 
nite pulses, in both time and frequency formulation, and 
use it to comment on general properties of the RABITT 
spectroscopy in the non-resonant case. 

In dipole approximation, the total hamiltonian 'H(t) of 
the target atom or molecule in interaction with a light 
pulse is given by a field-free component H plus an inter¬ 
action term [35], 

n{t) = H + F{t)o, o = i-d, (I) 

where F{t) = F{t)i is the external transverse light field, 
which for simplicity we assume to have constant polar¬ 
ization e, and O is an appropriate dipole operator. For 
example, in velocity gauge F{t) is the vector potential 
A{t) and O is proportional to the total canonical elec¬ 
tron momentum P = X]i=i Pij O = aP, with a being the 
fine-structure constant, while in length gauge F{t) is the 
light electric field E{t) = —adtA(t), and O is minus the 
dipole moment of the system, O = —jl. The wavefunc- 
tion IV’(t)) for the system, initially in the ground state 
\g), Ho\g) = \g)Eg, is 

\m) = \9)-^ f dt'F{t')Oj{t')m')), (2) 

where 0/(t) = exp{iHt)0 exp{—iHt) is the dipole oper¬ 

ator in the interaction picture. Unless stated otherwise, 
in the following we will use atomic units throughout. The 
r.h.s. of (2) can be expanded iteratively to arbitrary or¬ 
der in the interaction term [36], 

00 

= ( 3 ) 

n—0 

\i’^^\t)) = \g), (4) 

|^(n+i)(^)) = _,y dt'Eit')Oiit')\^^^\t')). (5) 
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The question now arises as to whether a truncated version 
of the expansion in (3) can adequately describe attosec- 
ond pump-probe experiments. The answer depends on 
the intensity of the laser, its duration and the coupling 
strength between all the states involved. With the strong 
VIS and IR dressing pulses available today, the contribu¬ 
tion of terms beyond lowest order may indeed become 
important [11, 37-43]. Rabi oscillations, for example, 
require the summation of the perturbative series to suffi¬ 
ciently high order to be reproduced across any given finite 
time interval. If both pump and probe ultrashort pulses 
are weak, however, the lowest-order approximation can 
be used to make accurate predictions. This is certainly 
the case of the RABITT technique described in the intro¬ 
duction, for which the probe intensity is kept small on 
purpose. The lowest perturbative transition amplitude 
from the initial ground state \g) to 
a final state ]/), Ho\f) = \f)Ef, featuring both pump 
and probe contributions, appears at second order, which 
can be cast in the following expression, 

Af^ =-i ff Fit 2 )F{ti) x 

J J — OO ( 6 ) 

X {f\0G+{t2-h)0\g), 

where we introduced the retarded Green function for 
the time-dependent Schrodinger equation of the field-free 
system, G+(t) = —i9{t)exTp{—iHt), with 9{x) being the 
Heaviside step function. We indicate the energy of a field- 
free state jj) indifferently as either Fi or Wi(= Fi/h), and 
energy differences as = coi — ujj. Equation (6) has a 
well known equivalent frequency counterpart, 

/ OO 

dujF{ufg - uj)F{uj)Mf^ (w), (7) 

-OO 


The FT of the individual pulses can be decomposed in 
the sum of an absorption (-|-) and an emission (—) com¬ 
ponent, 

E„M=F+M + E-(a;), (11) 

( 12 ) 

where the function (5 t(w), proportional to the FT of the 
characteristic function of the [—T/2,r/2] time interval, 

Sriuj) = [9{T/2 - |t|)] (a;) = (13) 

v27r TTW 

is a representation of the Dirac delta function. When 
replacing expressions (10-13) in (7), there are a limited 
number of contributions for which the frequency compo¬ 
nents from the two convoluted field FT overlap and which 
thus need to be considered, 

X (14) 

CX(T 0(t' 

/ OO 

du) Sriojfi -I- (j'uj0 — w) 5 t{oj -I- ctwq,) M fi{u}), 

-OO 

where cr = =Fl stands for photon absorption and emission, 
respectively. The last integral becomes negligible as soon 
as the energy-preserving condition is not satisfied, jw/i -I- 
auJa + cr'ujp I ^ 1 /T . If the two-photon matrix element 
is almost constant for |w -|- ctWq,] < 1/T, then 
we can replace it with the constant term Mfi(—au}a) 
and move it out of the integral. Using the convolution 
theorem, f f{x — uj)g{uj)duj = fg{x), together with 
Eq. (13), the remaining integral can be evaluated as 

J do; (5 t(A — w) (5 t(w) = (5t(A). (15) 


where F{u;) = {2 tt) f F(t) exp{iujt) dt is the Fourier 
Transform (FT) of the field, and (uj) is a two-photon 
transition matrix element, 

MfgH = (flOG+(u;g+u;)Olg}, ( 8 ) 

with the retarded resolvent G+(w) being defined as 

G+(u;) = J G+{t)F‘^^dt=(uj-H + iO+)-\ (9) 

a. The stationary regime. From equation (7), it is 
easy to derive the familiar formula for stationary transi¬ 
tion rates. To do so, let us suppose that the field com¬ 
prises a set of overlapping square pulses Fa(t), with dif¬ 
ferent frequencies uja and amplitudes Fh o, but all lin¬ 
early polarized along the z axis and having the same 
duration T, 


F(t)=^F„(t), 

06 

Fa{t) =Fa,Q COs{uJat + (fia) 9(T/2 - jtj). 


( 10 ) 


The transition amplitudes, therefore, becomes 

Ap} - J E ^“.0 E x 

Ct(T 0(7' 

X Sriujfi + crwa -I- (j'ujp), (16) 

Eq. (16) is the familiar stationary formula expressing 
the transition amplitude as a coherent sum of contribu¬ 
tions from individual time-ordered Feynman diagrams. 
When Eq. (16) is valid, it is possible to define a transi¬ 
tion rate W/i = limT^oo \Ap‘p^/T on account of being 
limT^oo 2-k 5'^{uj)/T 5{uj). 

Equation (16) differ from (6) and (7) in a fundamental 
way. While either expression (6) or (7) are applicable 
in the presence of intermediate resonant states, the sta¬ 
tionary expression (16) generally is not. The reason is 
that, the closer one gets to the resonance, the longer the 
exposure time required to legitimately factor out the two- 
photon matrix element from the integral in (14). Thus, 
for pulses comparable to or shorter than the characteris¬ 
tic lifetime of the resonance, equation (16) is not applica¬ 
ble as such, even if the truncated perturbative expression 
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is valid. In this latter case, despite the transition being 
second order, a stationary regime is never achieved and a 
transition rate cannot consequently even be defined. Fur¬ 
thermore, for long exposures the second-order transition 
amplitude may become so large that higher-order terms, 
possibly infinitely many of them, are required to achieve 
a physically meaningful result. 

b. The pump-probe scheme. Let us now examine 
how a finite-pulse formulation of the second-order transi¬ 
tion amplitude, such as Eq. (7), can be used to describe a 
pump-probe process. In pump-probe experiments the to¬ 
tal external field is the sum of a pump field Fi(<), which 
can be assumed not to depend on the pump-probe time 
delay, thus defining the time reference, and of a probe 
field F 2 {t',T) = F 2 {t — t) delayed with respect to the 
pump by a time lapse r (see Fig. la), 

F{t)=F^{t)+F2{t-T). (17) 

The FT of the total field has a simple parametrization 


E(t) 


(a) 


‘|E(a))| 

i 

(b) 
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'arg E((b) 
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the pump, 2 and 2 for the probe), replacement of (18) 
in (7) gives rise to sixteen terms associated to all possi¬ 
ble time-ordered two-photon transitions: 21, absorption 
of a pump photon followed by the absorption of a probe 
photon; 21, absorption of a pump photon followed by the 
emission of a probe photon; 12, absorption of a probe 
photon followed by the absorption of a pump photon, 
and so on. For example, the total transition amplitude 
for the absorption of one pump and one probe photon 
comprises two terms, 

(19) 

nOC 

■^fg = “* / duj F2{ujfg - uj;T)Fi{uj)Mfg{u}), ( 20 ) 

Jo 

poo 

^fg = “* / duj Fiiujfg - w)F2{uj;T)Mfgiuj), ( 21 ) 

which correspond to the time-ordered diagrams where 
the pump photon is absorbed first and last, respectively. 
Let us consider the first case in more detail. We can ex¬ 
pand the resolvent G^{ujg -bw) in the two-photon matrix 
element in terms of the generalized eigenstates |'0ae) of 
the field-free system, FI\ipas) = \i^ae)s, where a is a col¬ 
lective set of quantum numbers, other than the energy, 
sufficient to identify any given state (channel index). 




G f^aeGae.g 

UJg -\- U) - € 


( 22 ) 


FIG. 1. (Color online) Pump-probe scheme, (a) Temporal 
perspective: the second pulse is centered at a time r from 
the first pulse, which defines the time origin, (b) With non¬ 
overlapping pump and probe spectra, the photon distribution 
of the pump-probe sequence does not depend on the delay be¬ 
tween the two pulses, (c) The relative phase between different 
frequency components of the field, however, does depend on 
the time delay, (d) Since the same final energy can be reached 
with different combinations of the energies contained in the 
two pulses, the corresponding amplitudes can interfere con¬ 
structively or destructively depending on their mutual phases 
and, in turn, on the pump-probe time delay. 

in terms of the FT of the pump pulse and of the probe 
pulses at zero time delay (Fig. lb,c), 

F{oj) = Fi{oj) + F2{oj)e^‘^^. (18) 

In a two-photon transition with finite pulses, the energy 
preserving condition u!i-\-uj 2 = uifg is satisfied by several 
different pairs of frequency components (a;i,a; 2 ), which 
result in separate contributions that interfere to give rise 
to the total transition amplitude (Fig. Id). Changing the 
time delay between pump and probe pulses alters the rel¬ 
ative phase between all these different contributions, thus 
affecting the total amplitude, which becomes a function 
of r. 

If we consider separately the positive and negative fre¬ 
quency components of the field (1 and 1, respectively, for 


If I/) is either a discrete state or a generalised state be¬ 
longing to a featureless continuum (far from thresholds 
and from any resonant state), and the intermediate states 
contributing to (22) are either similarly featureless con- 
tinua or discrete states far from the resonance condition 
(virtual excitations), then A4fg(oj) is a smooth function 
of u! and, for sharply peaked field spectra, one recovers 
the familiar quasi-stationary expression for Aj:g as a fi¬ 
nite sum of weighted Feynman diagrams. In the presence 
of intermediate resonant states with lifetime comparable 
to or longer than the duration of the light pulses involved, 
however, Aifg(uj) has a sharp dependence on uj and the 
transition never achieves a stationary regime. In this 
latter case, the folding with the field in Eq.(7) must be 
evaluated to the full. 

c. The on-shell approximation. It is worth examin¬ 
ing the special case for Eq. (22) in which both the in¬ 
termediate states \ae) and the final state |/) = \(iE) 
are elastic-scattering featureless continuum states cor¬ 
responding to a same parent ion. In this case, the 
largest contribution to the two-photon transition ampli¬ 
tude comes from the intermediate states that are degen¬ 
erate or almost degenerate with the final state. This cir¬ 
cumstance is evident if the continuum states are approxi¬ 
mated with plane waves, which is a common assumption 
for energetic photoelectrons in multiphoton transitions 
(this approximation is employed in disguise, for example, 
in the strong-held [44-46] and in the soft-photon [47, 48] 
approximations, both of which are known to work well 
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sufficiently above threshold). Indeed, since plane waves 
are eigenstates of the dipole operator in velocity gauge, 
the only non-vanishing dipole transition matrix element 
is the one between two identical plane waves. 


{k\f\k') = kS’-^Hk-k') = k5^^\k-k') ^^^ ( 23 ) 

v2E 

Notice that such approximation applies when estimating 
the two-photon transition matrix element from a bound 
state \g) to the continuum, 

{k\p,G*iE, + ^)vAg) c + 

It does not imply, however, any net absorption or emis¬ 
sion of photons between free-electron states, which is and 
remains a prohibited process. We will call on-shell ap¬ 
proximation the assumption that the transition matrix 
element between unstructured continuum states is diag¬ 
onal in the energy, 

{PE’\0\aE) ~ Opa.(E)S(E - E'), (25) 

where Oap(E) is the integral of the actual transition am¬ 
plitude OpE,ae in a narrow energy interval e € (E — S, E-\- 
S) to which significant values of Oap(E) are hopefully re¬ 
stricted, 

pE-\-S 

Opo.{E)= / {l3E\0\ae)de. (26) 

Je-5 

The on-shell approximation is quite acceptable even 
when considering radiative transitions between the 
Coulomb or shifted Coulomb waves commonly encoun¬ 
tered in atomic ionization, and it becomes increasingly 
more accurate as the electron energy increases. For 
example. Fig. 2 (see also [49]) shows the continuum- 
continuum transition matrix elements in the hydrogen 
atom from three selected initial scattering states with 
£ = 0 and energies Eg =1, 2, 3 a.u., to £ = 1 scattering 
states as a function of the energy Ep of the final states. 
It is clear that the transition amplitudes are strongly 
peaked at Ep = Eg. In conclusion, using the on-shell 
approximation, the non-resonant two-photon transition 
matrix element from an initial state \g) to a final contin¬ 
uum state \I3E) through intermediate continuum states 
jaer), = {fdE\OG+{Eg -\-u})QaO\g), where 

is the projector on the intermediate continuum a, can be 
written as 


M 


(a) 

PE,g 


(w) 


Opa(E)OaE,g 
Eg Lj - E i0+ ' 


(27) 


The two-photon transition matrix element M.fg(oj) has 
thus assumed the form of a rational function which, apart 
for a simple pole in the lower half of the complex plane, 
depends only weakly on the frequency w. In the next sec¬ 
tion we will see that, with some additional approxima¬ 
tions, A4fg(uj) can be cast in a form similar to Eq. (27) 



FIG. 2. Absolute value, in atomic units, of the ex¬ 
act analytical reduced velocity-gauge dipole matrix element 
('i/'Bp ||Oi Il'i/’Bs) of the hydrogen atom, from three selected s 
scattering states (Eg = 1, 2, 3 a.u.), to several p states in 
the continuum. The sharp localisation of the amplitude at 
Ep ~ Eg underpins the validity of the on-shell approxima¬ 
tion. For more details, see [49]. 


even in the presence of intermediate and final resonance 
states. When this is the case, folding with the field com¬ 
ponents, as in Eq. (20), can be computed analytically for 
certain shapes of the light pulses. In the following, we 
will examine the relevant case of Gaussian pulses. We will 
subsequently apply the formula to the case of the non¬ 
resonant RABBIT transition and examine the effect of 
finite pulse duration on the RABITT beating frequency. 
The more general case of intermediate and final resonant 
states will be treated in Sec. III. 

d. Gaussian pulses. The vector potential of an ul¬ 
trashort laser pulse can be conveniently approximated 
with a linearly polarised Gaussian pulse A(t) = zA(t) 
parametrized as 

A{t) = cos[ajo(£ - to) + ^], (28) 

where Aq, wq, to, u and (j) are the amplitude, carrier an¬ 
gular frequency, central time, spectral width and carrier- 
envelope phase of the pulse, respectively. Several Gaus¬ 
sian pulses can be combined to give rise to arbitrary 
pulse sequences, or to chirped pulses. The absorption 
and emission components in the FT of a single Gaussian 
pulse, A(u!) = A+(a;) -I- A~(uj), are 

A^(uj) = — . (29) 

2tT 

As shown in App. A, in the case of Gaussian pulses, the 
folding of the field with a simple-pole function (u} — Zo)~^, 
as in Eq. (20), can be expressed in closed form. 


EjtklpE.g - U};t)Fi(uj) 
-oo Eg -f w — E -\- iO"*" 


iE^\T)E‘^^^w(zf), 

(30) 
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where is a form factor of the pulse sequence 




4cricr2 


X exp 


1/(52 r2 <J 2 S T 

— — I —2“*”^*- 

o-f (7i a at 


(31) 


with a = yja\ + cr|, cr* = Y CTj^ ^ + (72 and (5 = i5g + 

wi + a ;2 — E, while the complex parameter z|l is defined 
as 


,21 _ 

^ V2 




— E + Eg , 


(32) 


with w{z) = e~^ erfc(—iz) being the Faddeeva special 
function. The transition amplitude therefore, 

takes on the form 


^fE,g = ^''(T)e“=" ^a;3„(F;)0„£;..z«(4')- (33) 

Ol 


This last equation is essentially equivalent to the one for¬ 
mulated by Ishikawa and Ueda in terms of the Dawson 
integral (compare with Eq. 2 in [50]). 

In the region where the pulses do not overlap, the two- 
photon transition amplitude vanishes. How it gradually 
decays as a function of the pump-probe time delay is dic¬ 
tated by the product E^^{t)w{z^), which falls off like a 
Gaussian for jrj at- The left panels of Figure 3 illus¬ 



FIG. 3. (Color online) Spectrum for the two-photon ioniza¬ 
tion of the hydrogen atom from the ground state by means of 
a single (left panels) or a train of (right panels) XUV Gaussian 
pulses, in association with a 10 fs long IR probe pulse. Top 
panels (a,b): comparison between the energy integrated signal 
of the upper sideband, as a function the pump-probe time de¬ 
lay, computed ab initio (black dotted line) or with the model 
(blue solid line). Middle panels (c,d): energy and time-delay 
resolved spectra from the ab initio calculation (only states of 
even symmetry are shown). Bottom panels (e,f): energy and 
time-delay resolved spectra computed with the model. 


trate the photoelectron spectrum of the hydrogen atom 
ionized from the ground state with a single XUV Gaus¬ 
sian pulse with duration of 5 fs and central frequency 
40.8 eV, in association with an 760 nm IR probe pulse 
10 fs long, with an intensity of 10GW/cm2, as a function 
of the pump-probe time delay. The spectrum in the cen¬ 
tral panel is obtained ab initio by solving the TDSE for 
the atom in a numerical basis, while the bottom panel is 
computed using Eq. (33). The spectrum computed with 
the model, which includes all the terms proportional to 
the intensity of the probe laser, accurately reproduces all 
the features in the real energy-integrated (Fig. 3a) and 
energy-resolved (Fig. 3c) spectrum. 


H2n+r 
S^n “ 
H2n-r 






At 
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FIG. 4. (Color online) Quantum paths contributing to a side¬ 
band signal in RABITT spectroscopy. The amplitudes of both 
paths 1 and 2, in which one IR photon is absorbed, are mod¬ 
ulated by a phase factor , while those of paths 3 and 

4 are modulated by a phase factor As a result of 

the interference between the four amplitudes, therefore, the 
sideband signal beats with angular frequency 2ajiR. 


e. Red shift of RABITT beating with finite pulses. 
In RABITT spectroscopy, the amplitude of each sideband 
SB 2 n is given by the sum of four time-ordered two-photon 
amplitudes, 

A _ 2|H2n-lIR I 2|IRH2n-l I 2|H2„ + lIH. I 2|mH2„+l /0/|\ 

•^SB2„ — •^SB2n W "^SB2n •^SB2 ti •^SB2r, J 

where m indicates the emission component (negative fre¬ 
quency) of the IR pulse (see Fig. 4). In the limit of 
long pulses, the frequency of the RABITT beating is 2a;iR 
and the atomic phaseshift in the standard expression for 
the sideband intensity, /sb 2 „ = lo cos(2a;iRr -I- Acfa + 
A(/Jat) [15], is Atpat = arg[7W/i(a;2„-i) -b Mf^{uJtR)] - 
arg[AI/i(a; 2 n-i-i) + AI/i(— Wir)]. This is still the case if 
only one of either the APT or the IR has a long du¬ 
ration. Indeed, the RABITT frequency is given by the 
sum of the frequencies of the absorbed and of the emit¬ 
ted IR photons. If the XUV train comprises only multi¬ 
ples of the fundamental frequency Wir or if the probe 
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pulse is monochromatic with frequency Wir, then the 
only possible outcome for the RABITT beating is 2wia. 
On the other hand, when both the APT and the IR 


2(i)-25 2(0 2(0+5 2(0+25 



FIG. 5. (Color online) The XUV + IR above-threshold ion¬ 
ization amplitude is inversely proportional to the frequency of 
the IR photon. When both the XUV train and the IR pulse 
have finite duration, therefore, the signal is biased in favor of 
the low IR-frequency components. As a result, the spectrum 
of sideband beating in RABITT is red-shifted compared to 
the nominal 2LUin value. 


have finite duration, the RABITT beating is red-shifted 
with respect to the nominal 2a;iR value. This is be¬ 
cause non-resonant two-photon matrix elements, which 
are dominated by contributions from virtual states at 
the same energy as the final state, are inversely propor¬ 
tional to the energy of the last-exchanged IR photon, 
ME,i oc {Ei+uj 2 n±i — E 2 n)~^ = ■ Therefore, of the 

many IR wavelengths that contribute to the transition 
with finite pulses, long ones weigh more, thus biasing the 
RABITT beating towards the red (see Fig. 5). The right 
central and bottom panels of Figure 3 show the compari¬ 
son between ab initio and model calculations in the case 
of the RABITT ionization of the hydrogen atom, where 
a 5 fs long (fwhm) Gaussian APT formed by Gaussian 
XUV pulses with central frequency of 40.8 eV and dura¬ 
tion of 250 as is used in association with a 760 nm, 10 fs, 
lOGW/cm^ probe pulse. For these pulse parameters, the 
Fourier transform of the sideband oscillation in Figure 3b, 
reveals a beating frequency which is red-shifted with re¬ 
spect to the nominal RABITT frequency by an amount 
of 0.021 eV for the ab initio, in good agreement with the 
value of 0.017 eV predicted by the model. Part of the dif¬ 
ference between these two values is explained by the use, 
in the ab initio calculation, of a probe pulse with cosine- 
squared instead of Gaussian envelope, which permits us 
to reduce the size of the quantization box. 


III. TWO-PHOTON RESONANT MODEL 

In this section we will use Fano formalism to compute 
two-photon transition amplitudes for the case in which 
the intermediate and/or final continuum states feature 
isolated resonances. After a short overview of Fano’s 
main results and a discussion of the phase properties of 


the one-photon Fano transition matrix element, which 
will be relevant for the following of this section, we will 
consider two-photon transition amplitudes for the case of 
single-channel continuum states with at most one isolated 
resonance. The generalisation to multiple single-channel 
isolated resonances will be straightforward. At the end 
of this section, we will discuss the extension of the model 
to the multichannel resonant case. 


A. Phase properties of Fano transition amplitude 

In stationary conditions, metastable states manifest 
themselves in single photoionization spectra as charac¬ 
teristically asymmetric peaks [51]. The asymmetry of 
experimental resonant profiles can be explained with 
the well-known Fano’s formalism [52]. In the simplest 
formulation of Fano’s approach, the field-free hamilto- 
nian H is given by the sum of an unperturbed compo¬ 
nent i7o and a “configuration interaction” component 
V, H = Hq + V, where the eigenstates of Hq com¬ 
prise a featureless continuum je) and a bound state ja), 
Ho\e) = ]£)£, Ho\a) = \a)Ea, while the configuration in¬ 
teraction only couples the bound state to the continuum, 
Vae = {a\H — Ho\e). If the coupling I4e depends only 
weakly on the continuum index e, the continuum eigen¬ 
states of the full hamiltonian, HIiIje) = \4’e)E, can be 
expressed as 


I'fpE) = \E) + 


^) + J de\e) 


14 , 


E — e + i0+ 


VaE 

E-E 


■, (35) 


where Ea{E) = Ea+l2\a{E)—iTa{E)/2 is a complex func¬ 
tion of the energy, with the so-called energy shift Aa{E) 
and width Ta{E) defined as Aa{E) = P^delVeal"^/{E — 
e) and Fa(i?) = 27r^jt4;aP- The pole of E — Ea{E) in the 
negative complex plane is, by convention, the complex 
resonance energy. Notice that the solution (35), which is 
readily obtained by projecting the Lippmann-Schwinger 
equation 1'0£;) = \E) + [E)V\'iIje) on the basis of un¬ 
perturbed states, differs from Fano’s original solution by 
a complex normalisation factor. Here we will use the 
form (35) because it is normalised, {i1^e\'4’e') = S{E—E'), 
and continuous with respect to E. The energy shift 
Aa{E) and width ro(A) depend only weakly on the en¬ 
ergy E, so that one can assume they are constant in 
the energy region of interest, Fa(A) ~ Fa(Ao) = Fa, 
Aa{E) ~ Aa{Ea) = A^. In these conditions, which 
we assume to be fulfilled, the complex resonance energy 
is thus well approximated as Ea — Ea + Aa — iFa/2. 
For our convenience, we will indicate the real part of 
the resonance energy as Ea = Re[ii'o] = Ea + Aa- It 
is customary to define a reduced energy variable e = 
2{E — Ea)/Ta, and a distorted continuum component 
ja) = ]a) -I- P '^d£\e)Ve,a/{E — e) which incorporates the 
original bound state ja). Using this notation, Eq. (35) 



























can be reformulated as 


\i’E) 


\E)^. + \h) 


e +1 


1 1 

ttVeo e + i 


Notice that for e —>• ±oo, |'0 _e) converges to \E). 
define a resonant phaseshift (j)E as 


(j)E = 'e/2 arctane, 



(37) 



which is a continuous, monotonically increasing function 
of E, with (/)_oo = 0, (/foo = 'E- If the \E) channel func¬ 
tions are used as reference asymptotes, 4>e is associated 
to the on-shell scattering matrix s{E) for the collisional 
excitation of the resonance, 


FIG. 6. (Color online) Trajectory in the complex plane of the 
77(e) resonant factor in Fano’s dipole transition amplitude, as 
the reduced detuning e varies from large negative to large pos¬ 
itive values. The three panels correspond to different values 
of the q parameter. See text for more details. 


s{E) = e^^^ = {e-i)/{e + i). (38) 

Conversely, e = — cos^/sin^ (we will drop the energy 
suffix from (j)E, when E is clear from the context). In this 
formalism, the dipole transition matrix element between 
an initial ground state \g) and a final resonant continuum 
\iPe) can be written as 

{^E\0\g) = OEg'-^^, Oes = {E\0\g), (39) 

where qag is a real parameter that measures the relative 
strength of the transition from the ground state to the 
autoionizing state, relative to that of the direct-ionization 
process. 


qag — 


O-aa 

TrVaEOEg 


(40) 


The resonant factor 77(e) = (e -I- qag)/{i + i) in the 
complex conjugate of the dipole transition amplitude (39) 
can be written as the sum of a constant term plus a 
second term proportional to the unimodular function (e— 
i)/{e + i), 


77(e) 


e + g 

e -I- i 



1 -\- iq e — i 
2 e -t- ^ 


(41) 


This means that, as e increases from —oo to -boo, 77(e) 
describes counterclockwise a circle in the complex plane. 
The circle is centered at (1 — fg)/2, it has radius r = 
\/l -b / 2, and it both starts and ends at 1, intercepting 
the origin at e = —q. Figure 6 illustrates the trajectory 
of i?(e), from large negative (e <C — 1) to large positive 
detunings (e ^ 1), for three representative values of q: 0, 
0.5, and 1. This geometrical interpretation can be visu¬ 
alized even more clearly by defining the angular variable 
ip = arctan(g) € (—7r/2,7r/2), equivalent to the one in¬ 
troduced by Ott et al. [53] in the context of the dipolar 
response of a Fano resonance, with which 77(e) becomes 

77(e) = re-*^ + (42) 


Indeed, the latter experiences a discontinuous jump of tt 
in correspondence of e = —q, when the circle intercepts 
the origin. 


arg77(e) = arg {e*‘^ | 


= 4> + arg [cos(^ -b ip)] = 
= arctan(e) — ■n0{e -b q). 


(43) 

(44) 

(45) 


Far from reflecting a real discontinuity in the physical 
properties of the system, this circumstance simply re¬ 
flects the fact that arg( 2 ;) is discontinuous at the origin. 

The photoionization cross section <j{E), which is pro¬ 
portional to the square module of the dipole transi¬ 
tion amplitude to the resonant continuum, therefore, is 
the product between a background smooth cross section 
o'bg(F'), associated to direct photoionization, and the cel¬ 
ebrated Fano profile, 


a{E)=a^^{E)^^±^. (46) 

In this latter expression, any information on the rela¬ 
tive phase between different frequency components of the 
photoelectron wavepacket generated by the interaction 
of the impinging ionizing light is lost. Therefore, while 
Equation (46) is sufficient to interpret one-photon ioniza¬ 
tion experiments such as those conducted at synchrotron 
facilities, when two or more photons are exchanged in 
a coherent transition, the relative phase of intermediate 
transition amplitudes becomes essential and we must go 
back to Eq. (39) instead. 


B. Resonant two-photon transition matrix element. 

To derive the analytical formula for finite-pulse res¬ 
onant two-photon transition amplitudes, we first need 
to obtain an approximated analytical expression for the 
two-photon ionization matrix element M.pE,g{'^)-, 


Notice that even if the resonant dipole transition ampli¬ 
tude (39) is a continuous function of e, its phase is not. 


M/3. 


bff(w) = ^de 


{lppE\0\'^aE) {'ipae\0\9) 
Eg -\- u) — S' 


(47) 
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To do so, we assume that the continuum branches in both 
the intermediate states, I'i/’as), and final states, \iI)pe), 
can be expressed, using Fano’s formalism, in terms of 
known bound and continuum eigenstates of a reference 

hamiltonian i?o, ^^o|a) = Ea\a), Hq]})) = -ffo|7£^) = 

e|7e), 


\i>aE) 


\aE) + 

m) + 



de\ae) Vge.a ^ 
E — e + i0+ / 



de\j3e) Vpefi \ 
E — £ + iO+y 




E-Ea 


,(48) 


Vb,0E 

E-Eb 


(49) 


where V denotes the field-free electron-electron interac¬ 
tion not included in Eq, e.g., 14,ae = (n|^ ~ ^^olae). 
The interacting-continuum wavefunctions in Eqns. (48) 
and (49) are normalized as (V’b'IV’-e) = ~ ^)- The 

suffixes a and (3 identify the ionization channel in the 
intermediate and final states, respectively, i.e., the col¬ 
lection of discrete quantum numbers needed to specify 
the asymptotic state of the parent ion, of the photoelec¬ 
tron, as well as their angular and spin coupling. Let us 
separate in AipE,g{^) the contribution of the interme¬ 
diate bound states {|n)}, g{uj), from that of the 

intermediate continuum states IV’as), 


expression 


M 


ib) 

PE,g 


(w) 


%no OpE ,no Ong,g 
^Eb ^ ^nog 


(56) 


The latter expression is applicable even in the case of 
multiple intermediate states that contribute to the transi¬ 
tion amplitude by means of virtual excitations and which 
are clustered in an energy region that is small if compared 
with the detuning w — uJnog from the absorption of the 
first photon. Let us now consider the contribution from 
the intermediate continuum states. Replacing eq. (53) 
in (52) we find 


M 


(c) 

PE,g 


(uj) 


^ea Qag ^ 

Eg -\- UJ — 6 -\- ^ 


(57) 


To advance further, we must find an expression for 
the continuum-continuum resonant transition amplitude 
{4’l3E\0\4’ae) in terms of a limited number of almost- 
constant parameters. In analogy with the Fano formula 
for the dipole transition from bound states, we first take 
out from this matrix element the term that involves only 
transition matrix elements between states in the unper¬ 
turbed continuum. 


Mi 


((c) 


(b) f, .\ i'^PE\0\n) Ong 


I / x _ \7/3J 


bJng + iO'b 


M 


(c) 

PE 


= / 


de 


{'lppE\0\'tl^ae) {lpae\0\g) 

Eq UJ — S' H“ 


(50) 

(51) 

(52) 


The transition matrix elements between a localised state 
and a Fano continuum can be accurately parametrised 
with Fano’s formula. 


i^Mg) = + (53) 

Cea - * 

{iPpEpln) = ^-^^±^OpE,n- (54) 

^Eb — * 

The contribution from bound intermediate states, there¬ 
fore, is readily written as 


J^?E,gU) = E 

n 


^Eb T Qjjji E^pEjTiE^ng 
eEb-i W - ujng -b i0+ 


(55) 


{tppE\0\tpae) = {tppElOllljae) — „ ^^'1. Oba 


E- E* £- Ea 


■ -b 


E- Et £- E„ 


= ('0/3£|O|V’ae) + - 


1 


Oba 


E £Eb ~ i Eea + * - 
OpE, 


EVb^pEVa 


+ , (fea + gab) + [cEb + qia) 

Vh.RE Vrvfn. 


(58) 


where the barred states represent the Fano continuum 
without the bound component, i.e.. 


IV'/sb) = \PE) 


del/de) Vpe,b Vb,pE 

E-£ + i0+ E-Eb' 


(59) 


and we used the relation F^ = 27r|I4,o!£;P- By applying 
the on-shell approximation, and assuming that Oap = 
Oap{E), Va,aE, and Vb^pE are sufficiently slowly varying 
functions of E, it is easy to show that 


In practical cases, this expression can often be restricted 
to the contribution from a limited set of intermediate 
bound states, or even from just one of them. For ex¬ 
ample, in the excitation of helium from the ground 
state to the doubly excited states with iV = 2, the biggest 
role in Eq. (55) is played by the intermediate ls2p state, 
for which the oscillator strength with the N = 2 states 
is very large and the background ionization amplitude 
is very small (gj,„ » 1. The 2p^ ^ ls2p ^ Is^ is a 
characteristic example). In this case, if the intermediate 
bound state is non-resonant, one can use the simplified 


iijpEm^e) = Op^6{E - e) + (60) 

^ 1 El Pol eEb ^ea 

TT £ - £) -b f0+ (Cea + i){eEb “ *) 

Indeed, to compute the transition matrix element be¬ 
tween the two modified continua, it is sufficient to close 
the integration path with a semi-circular path in either 
the upper or the lower half complex plane, where the 
argument of the integral decreases quadratically with 
respect to the integration variable, and apply Cauchy 
residual theorem. By combining Eqs. (58) and (60), the 
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dipole transition amplitudes between the two Fano reso¬ 
nant continua can be approximated as 

{ll^liE\0\'li)ae) = - e)-|- 

^ ^I3a ^Eb Cea ^ 

£ - E + i0+ n{esa + i){eEb - i) (61) 

Ob,ae{^ea-\-Qab) , O fS E ,a{^Eb-\-Qha) _ Obn. 

Vb,l3E yae,a TT Vae,a V^b,/3E 

T^ieea + i)ieEb - i) 

We can now insert this expression in the continuum con¬ 
tribution (57) to the two-photon matrix element. The 
integral of the argument proportional to a Dirac delta 
function is evaluated immediately, while the other two 
terms require a more careful discussion, 


M 


ic) 

l3E,g 


(uj) = 


^^a^aE,g ^Ea Qag 

Eq OJ — E ^Ea — ^ 


+/2+/3, (62) 


where 


1/^ f £ea T Qag £ea £Eb Ej^aOcte 

£Eb — i J £^ + 1^“^ + *0''’ £ — Eg — UJ — i0+ 


and 


h = 


/ 


^ae^gds 


Qag 


Eq -\- uj — s -\- 


Ob,ae{^ea~^Qab) 

Vb.dE 


eea - 'I 
0/3E,a{^Eb + Qhn. 

Vcce.a 


(63) 


Ob 


brVb.^E Vce.a 


7r(e£a + i)ieEb - i) 


For large values of e, the argument of the integral in E 
is inversely proportional to £^. Therefore, this integral 
can be conveniently computed by closing the integration 
circuit with a semi-circular arc in the lower half of the 
complex plane, provided that the transition matrix ele¬ 
ments are only weakly varying on the additional arc, for 
large enough arc radii. The result is 


h 


£Ea £Eb £Ea T Qag j^aEaE,g 

i^Ea ~ 0(ei5b ~ i) £Ea + * UJ — E -\- Eg 

£Eb T ^ Qag ^ E^c^Eoi^E,g 

£Eb — i £Ea + * UJ — UJag 


(64) 


The last integral, I 3 , has only one simple pole in the 
lower complex plane and hence it also would be con¬ 
veniently computed by closing the integration circuit 
in the lower half of the complex plane with the arc 
TR = {Re-^^,^&[0,7r]}, 


h 


lim 

R—¥0Q 


I I 3 (z)dz 


/ l3{z)dz 

'Tr 


(65) 


where X 3 (e) indicates the argument of the integral in (63). 
In contrast to the previous case, however, the absolute 
value of 13 (e) decays only as |£|“^. Instead of vanishing 
as i? —>■ 00 , therefore, the contribution of the arc con¬ 
verges to a finite value that must be taken into account, 
and which is easily computed (as usual, we assume that 


all the matrix elements are constant in a region of the 
complex plane sufficiently large to attain reasonable con¬ 
vergence of the circuit integral), 


lim / l 3 (z)dz = 

R-yco 


^E If ^O-eE olE ,g 

Vb^pE(£Eb — i) 


( 66 ) 


The value of the total integral I 3 , then, becomes 


h = 


X 


Eba Ea,g 

-I - 

14 ,/3 £Eb — i 

Eb^a (Qab 0 

Vb,/3 


ba Qag i Ea^g 
‘2 UJ UJbg £Eb f 


(67) 


E^^a(£Eb E Qbo) Eba/'^ 

bet,a ^b,0^a,a 


In summary, the expression for the intermediate- 
continuum contribution to the two-photon resonant tran¬ 
sition matrix element is 


M 


(c) 

dE,g 


(uj) 


£Ea T Qag £Eb i E0ciEct^g 

£Ea + * £Eb ~ i Eg -\- iu — E lO"*" 
Qag ^ £Eh “f ^ E^^EQ>^g ^ 

£Ea + i £Eb — i UJ — UJbg 


+ (Qab - i) 


Qag f ^a/^ Eb^aEct^g ^ 

£Eb — i Vb^pE UJ — UJ~ag 

Qag ^ EfH^aEa^g 


( 68 ) 


'^Va,a.(€Eb + Qla) 

eEb-l UJ - UJ'ag 

Qag ^ 14 ,ct EbaEa^g iEb^aEa^g 


£Eb — i Vb,p UJ — ujbg Vb,p(eEb — i) 


This approximate algebraic expression for the two- 
photon transition matrix element in the presence of both 
an intermediate and a final autoionizing states is one of 
the main results of the present work. 

/. Case of no final resonances. In the relevant spe¬ 
cial case in which there are no final resonances, Eq. (68) 
simplifies considerably since one can take its limit for 
vanishing radiative and non-radiative couplings involv¬ 
ing the \b) state. The result is 


M 


C) 

dE,g 


(uj) 


£Ea F Qag EjS^Ea^g 

£Ea F i Eg -\- UJ — E -\- i0“*" 


(69) 


/3a- 


1 


£Ea F * 


: (Q~ag - i) 


EjSaEo 


■^ag 


where the parameter j5a = t^E^^ aVaalEpa is a pure 
number that depends solely on the properties of the 
atomic system. When considering resonant two-photon 
transitions with long overlapping pulses with frequen¬ 
cies uJi and UJ 2 and duration larger than the lifetime 
of the intermediate resonance, Eq. (69) can be simpli¬ 
fied further, since the conservation of energies applies, 
E = Eg + uji + UJ 2 . With few algebraic passages, it is 
easy to show that the matrix element appropriate for the 
time-ordered diagram in which photon uji is absorbed 
first becomes 


M 


(c,21) 

dE,g 


(uJi) = 


EjSofEot^g CEia F Qag(^ Ta2) F i^a2 
UJ2 £Eia F i 


( 70 ) 
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where we introduced the new real parameter 


W2/3o _ (^l3,a 


(71) 


which measures the relative strength of two alternative 
paths for the dipole transition from the intermediate 
bound state |a) to the final continuum \PE): a direct 
one, Op^a, and an indirect one, OpaVaal ^2 hr which the 
transition is mediated by the non-radiative coupling of 
the bound state with the intermediate continuum \aE) 
followed by the dipole coupling between the intermediate 
and final continuum. Notice that in the formulation (70), 
the reduced energy term e is always relative to the en¬ 
ergy of the intermediate state reached from the ground 
state by the absorption of the first photon, exactly as in 
the one-photon formula (39). It is interesting, therefore, 
to analyze more in detail the similarities and differences 
between expression (70) and that for one-photon transi¬ 
tions. First of all, if we define an effective q parameter 
as 


giff ^ = gag(l - 7a2) + *7a2, (72) 

the resonant factor is formally the same in either expres¬ 
sions. 


I (21) 

'^(21) _ ^Eia ~r 9eff (73) 

^Eia + * 

Only in the case in which the intermediate bound state 
I a) is not radiatively coupled to the final continuum 
= 0 => = Qdg), however, do the resonant 

factors in the one-photon and the two-photon transition 
amplitude actually coincide in value and not in form only 
(see green line in Fig. 7). In general, if 7^2 ^ 0, is a 
complex number which depends on the frequency of the 
second exchanged photon. The resonant factor in the 
two-photon transition matrix element can also be writ¬ 
ten as 


7^(^l)=7,, + (l-7,,) ^^^° + g^^ (74) 

e^ia + ^ 

which is the same factor as in the one-photon case, scaled 
by (1 — 7 a 2 ) and shifted along the real axis by 7 q 2 . In 
particular, as the reduced detuning tEia is increased from 
—00 to - 1 - 00 , 72.i^^i(e_Eia) still describes counterclockwise 
a circle that starts and ends at 1. In contrast to the 
one-photon case, however, if 7^2 0 , the circle does not 

intersect the origin. In particular, if 7^2 < 0, the circle, 
which is expanded with respect to the one-photon case, 
intersects the real axis at 7^2 and at 1 , thus encircling 
the origin. This means that the phase of the two-photon 
transition matrix elements experiences a full 27r excur¬ 
sion. If, on the other hand, 7 a 2 > 0, the circle is con¬ 
tracted and it misses the origin. In this latter case, the 
phase of the two-photon transition matrix elements expe¬ 
riences a finite excursion but no overall variation (see blue 



-4 -2 0 2 4 


E 

FIG. 7. (Color online). Argument of the resonant factor 77^^ 
(73) of the two-photon matrix element (70) as a function of the 
reduced detuning of the pump photon from the intermediate 
resonance. Brown line: the ground state is radiatively cou¬ 
pled to the discrete but not to the homogeneous component 
of the intermediate state, i.e., q — >■ 00 . Green line: the ground 
state is radiatively coupled to both the discrete and homoge¬ 
neous components [q = 1), but the discrete component is not 
radiatively coupled to the final state, i.e., jSsa —t 0. Blue line: 
both the discrete and continuum intermediate components are 
radiatively coupled with the initial and final states. 


line in Fig. 7). Furthermore, since 7 a 2 is proportional to 
a; 2 , the sign of 7 q 2 for the emission of the second photon 
is the opposite of that for its absorption, the full 27r phase 
excursion and the no-net phase excursion case are both 
simultaneously present, one for the upper and one for the 
lower sidebands of the resonant two-photon transition. In 
the particular case in which 7^2 = 1 , the circular com¬ 
plex trajectory contracts to a point, = 1 , so that 

the two-photon amplitude does not bear any sign of the 
intermediate resonance (the amplitude for the opposite 
sideband, however, would still exhibit a pronounced 27r 
phase excursion). 

So far, we have considered only the case in which the 
photon close to the resonance is absorbed first. In fact, 
the same final state is also reached by the path in which 
the photon W 2 is exchanged first, and whose matrix el¬ 
ement is { 102 )■ In XUV-pump IR-probe experi¬ 

ment, where Wm <C Wxuvi nnd where the first excitation 
energy of the ground state is typically much larger than 
WiR, the contribution of the second path is generally small 
and it is often disregarded. Yet, the total transition ma¬ 
trix element should be computed as the sum of the two 
time-ordered contribution. If the path in which 102 is ex¬ 
changed first is not resonant, then we can imagine that 
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this term contributes with a small complex constant to 
the total transition. In principle, therefore, the inverted- 
order transition has an effect similar to that of 7 a 2 , as 
it shifts the transition amplitude trajectory towards or 
away from the origin. 

The limit in which only the intermediate state |a) is 
radiatively coupled to the ground, while the intermedi¬ 
ate continuum is not [qag = oo), is also interesting, since 
it effectively reproduces the assumptions that have been 
made in [29] and [31], and which lead to a neat tt ex¬ 
cursion of the transition amplitude phase, as shown by 
the brown line in Fig. 7. In the general case, where the 
discrete-continuum dipole coupling is not negligible, the 
typical abrupt tt discontinuity of (42) disappears (see 
blue line in Fig. 7). 

g. Case of no intermediate resonances. Two-photon 
excitation of a metastable state in the final continuum, 
with no intermediate resonances, which has been ex¬ 
plored in the past by Cormier et al. [54], is a second 
relevant case. The frequency-dependent two-photon ma¬ 
trix element for this case is readily obtained from the 
general formula (68) by suppressing all the terms that 
involve the intermediate state jo). 


M 


C) 

PE,g 


(w) 


eEb “b ^ CiC)t;)^(^(D 
CEb — iEg+uj — E + f0+ Vb^jsieEb — *) 


If we specialise this formula to the long-pulse limit, and 
assume the conservation of energy E = Eg + uji + uj 2 , we 
obtain 

= (75) 

^ Jo P — 7. 


where we introduced the new real parameter 


7b2 


W20&,a 

Vb^pOpa 


(76) 


A first surprising aspect of the resonant transition matrix 
element (75) in the present model is that it has a purely 
imaginary q parameter, g = z(l -|- 762 )■ As mentioned at 
the beginning of this section, however, when autoionizing 
final states are involved, the contribution of intermediate 
virtual bound states can be very large and, when added 
to (75), they give rise to an effective complex q parame¬ 
ter with comparable real and imaginary components, as 
predicted in [54]. Notice that, if the radiative coupling 
between the intermediate continuum and the final bound 
state is sufficiently large, it is in principle possible to se¬ 
lect a value of W 2 such that 1 + "/b 2 vanishes, thus making 
the transition amplitude disappear at one of the final res¬ 
onances, as it happens at e = —q for a traditional Fano 
profile. 


C. Time-resolved transition amplitudes 

From the expressions for the continuum (68) and dis¬ 
crete (55) contribution to the resonant two-photon tran¬ 
sition matrix element, we can now proceed to compute 


the full transition amplitude associated to a pair of Gaus¬ 
sian pump and probe pulses. To do so, we will fold the 
transition matrix element with the FT of the field, 

as prescribed in Eq. (7). Except for the last term in (68), 
which does not depend on the integration frequency vari¬ 
able w, all the other terms in either (68) or (55) depend 
on w through elementary factors of the form (oj — wo)~^- 
The folding in (7), therefore, can easily be carried out 
using Eq. (30). In the case of the absorption of photon 
1 followed by that of photon 2, the expression for the 
transition amplitude reads 


g = E^^Wpo.O^E,g X 

^ CEb-I 

X / ~^'^°g w(4i) -p {q^^g - z)w(z| ) X 


A 


eEa + I 

eEb + q-ba 


I 


+ 


eEb + i eEa + * 

21 

E (Tt Vb,pOpa{eEb + i) 


^ba{qab ^) Cba 
eEb + i 


E 


eEb + Qbn OpE.nOng 


21 




eEb + i OpaOaE,g I 

where we introduced the additional parameters 


(77) 


_ ra/2 Ob,a _ K.g Oba 

Vb^pE Opa ’ Vb,p Opa 


(78) 


Equation (77), which is one of the major results of this 
paper, depends on a minimal number of parameters for 
the radiative and non-radiative couplings between all the 
essential states involved in the dynamics, as well as the 
parameters of the pump and probe impinging pulses, in¬ 
cluding their time delay. Once the parameters of the 
model are established, therefore, this formula is able to 
provide, at a negligible computational cost, full energy 
and time-delay resolved attosecond pump-probe photo¬ 
electron spectra in the presence of both an intermediate 
and a final resonance for arbitrary pairs of (weak) pulses. 
Eurthermore, this result is trivially extended to the case 
of an arbitrary number of Gaussian pulses, to represent, 
e.g., the effect of an attosecond pulse train, as well as to 
an arbitrary number of isolated resonances either in the 
intermediate or in the final states. 

It is now interesting to consider more in detail the case 
of no final resonances, for which the transition ampli¬ 
tude (77) simplifies to 


^PE,g = ^ {r) C'^^^Opo.O^E.g X 




eEa 






O pE,nOng ( 21 \ 

^ 13 a. J^OiE^g 


( 79 ) 
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in relation to the simple two-photon transition matrix 
element (69) discussed earlier in this section. In partic¬ 
ular, we want to examine the effect of using finite pulses 
on the complex trajectory of the two-photon transition 
amplitude as a function of the central energy of the pump 
pulse. Each of the upper panels in Fig. 8 shows the com- 



FIG. 8 . (Color online) Complex trajectories of resonant 
finite-pulse two-photon two-color absorption amplitudes (up¬ 
per panels), and corresponding phase variation (lower pan¬ 
els) as the reduced intermediate energy detuning in¬ 

creases from large negative (—30) to large positive values 
(-1-30). In each of the upper panels, increasingly shorter pulses 
(utFa = oo, 2, 1, 0.5) give rise to progressively more con¬ 
tracted trajectories (traced with lighter color). Each column 
has a different set of q and 7 parameters (see text for details). 


plex trajectories of the transition amplitude (79), with 
^Eia S [—30,30] and at four different pulse durations, 
atTa = 00, 2, 1, 0.5, for a selected pair of resonance pa¬ 
rameters: a) gas = 20 » 1, 702 = 0, b) gag = 1, 7o2 = 0, 
c) gag = 1, 7a2 = 0.2, d) gag = 1, 7a2 = -0.2. The 
amplitudes are normalised so to start at the same ref¬ 
erence point on the real axis, which corresponds to the 
asymptotic background transition amplitude. The lower 
panels, Fig. 8.e-h, show the transition amplitude phase as 
a function of the reduced pump detuning e Ea ■ When the 
transition amplitude is dominated by the contribution of 
the intermediate autoionizing state (gag ^ 1, Figs. 8a, e), 
the observed phase excursion is always tt. The shorter 
the pulse duration, the wider the step (Fig. 8e). If g is 
finite but the intermediate autoionizing state is not ra- 
diatively coupled to the final continuum, the trajectory 
intercepts the origin, but only in the limit of long pulses, 
while, for short pulses, folding with the pulse spectra 
contracts the circular trajectory towards the asymptotic 
background value (Fig. 8b). In particular, the phase loses 
its discontinuity, giving rise to a sigmoidal profile with 
no net phase change, with features that are progressively 
less pronounced as shorter pulses are employed (Fig. 8f). 
The effect of finite pulses, therefore, is similar to that of 
a direct dipolar coupling between the bound state and 
the final continuum or, as we will see later in this sec¬ 
tion, to that of multiple intermediate channels. A similar 
dependence on pulse duration is observed for 7^2 > 0 


(Figs. 8c, g). The complementary case of 7a2 < 0 (if 
7a2 > 0 for probe absorption, 7^2 < 0 for probe emission, 
and vice versa) (Figs. 8c, g) is more interesting because, 
in the long-pulse limit, the phase experiences a full 27r 
jump, transitioning to the continuous excursion with no 
net phase through a point, for a definite finite pulse du¬ 
ration at, at which the phase has a discontinuous jump 
or, stated otherwise, at which the resonant two-photon 
transition amplitude exactly vanishes. 

If the energy of the second photon is much larger than 
the natural width of the intermediate resonance, a;2 ^ 
Tq, Eq. (79) can be further simplified. In fact, if the 
more stringent assumption |e/a| ^ go holds, we recover 
the expression given in Eq.(6) of [10], 




-i{u)2'r+rl>l+4>2) , 


w{Ze) + (/3a - iEi)i^ag - i)wizfj 


(80) 


which was indeed justified in the context of helium pho¬ 
toionization in the region of the doubly excited states 
converging to the N = 2 threshold. 

h. Correspondence between intermediate-energy scan 
and final-energy resolved photoelectron spectrum. So 
far, when commenting the case of no final resonances, 
we have concentrated our attention on the phase of the 
resonant two-photon transition amplitude as a function 
of the central frequency uji of the pump pulse scanning 
the resonance, for a given value of the final-energy detun¬ 
ing 5 from the nominal value Eg-\-uji-\-uj 2 . Alternatively, 
one can keep wi constant and study the dependence of 
the transition amplitude on the final energy instead. In 
either cases, the variation of the amplitude is essentially 
dictated by the resonant argument of the Fadeeva 
function in (79), 

g. - 

[uJi- LO-ag) + -\{E - Eg - LOi-LO 2 ) - i— , 

?81) 

all the other terms in (79) having, in comparison, only 
a weak dependence on E and wi. From Eq. (81), the 
similarity between these two cases is evident: in the first 
case (scan over wi), the second term in parenthesis is 
constant while the first increases linearly with wi; in the 
second case (scan over E), the first term in parenthesis 
is constant while the second increases linearly with E. 
In attosecond pump-probe experiments, furthermore, the 
pump pulse is oftentimes much shorter than the probe, 
and hence a\la'^ ~ 1. In these conditions, therefore, the 
two cases become essentially equivalent. 

i. Monochromatic limit. It is instructive to ascer¬ 
tain that the formula for the finite-pulse resonant two- 
photon transition amplitude (79) approaches the station¬ 
ary expression (70) in the limit of long overlapping pulses, 
i.e., assuming pulse durations much longer than the res¬ 
onance lifetime, atTa ^ 1, and time delays negligible if 
compared to the duration of the pulses, t at. For 
at 00, the argument of the Faddeeva function, z, 
tends to (-boo, 0) € C, so one can use the first term 


21 ^ 

J2 
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in the asymptotic expansion of w{z) restricted to the 
real axis, w{x) ~ a: £ IR [55]. For ~ 

ft /{Ei + wi — Saf — i?), the Faddeeva function 
becomes w(z|l) ~ i^j2/TT a^^{Ei +a;i — dalja'^ — E)~^. 

If we neglect the effect of intermediate bound states, we 
obtain 

jE-i „_ I _ egg + gag 

Eg+wi- dal/a"^ - E esa + i 

Eg + bji — 5a1/u'^ — Ea^°^ ^ 

At the nominal energy of the transition (<5 = 0), and 
using the energy-preserving condition E = Eg + uji + u} 2 , 
the transition amplitude finally becomes 


oc 


I (21) 
egia + ggff 

Cgia + i 


(82) 


elements involving the continua are smooth and slowly 
varying functions of the energy, however, generalisa¬ 
tion to the case of an arbitrary number of intermedi¬ 
ate and final continua is straightforward. As shown in 
Section 4 of the original Fano paper [52], the case of 
a bound state ja) coupled to several unperturbed con¬ 
tinua Jae), Vaa = (aejiJja), can be reduced to that of 
the bound state ja) coupled to a single “resonant” con¬ 
tinuum li?e), VRa = \/j2a \^aa\^ = {R£\H\a), plus a 
set of fully decoupled residual featureless continua Ja'e), 
(Q;'eji7ja) = 0, by means of a unitary transformation of 
the degenerate unperturbed continua. 


ii?£) = E \cX£^UaR^ jd s') — ^ ] \o^s)U(yct'^ 



as anticipated. 

j. Case of non-overlapping pulses In the presence of 
an intermediate resonance ja), instead of plummeting as 
soon as jrj > at, as it was the case for non-resonant tran¬ 
sitions (compare with Sec. II Od), the two-photon signal 
persists even for t > at, decaying exponentially as 
(tq = F”^). To see this, let us consider the transition 
amplitude at zero final energy detuning, when t ^ at, 
and the first photon absorption is right on resonance, 

wi = Re[wog], 


i f at T 
^\2Ta at 


(83) 


[/tf/ = f/[/t = 1, (87) 

Furthermore, the residual decoupled continua \a's) can 
be chosen so that only one of them, which we will call 
\D£), is radiatively coupled to the ground state, Oog = 
{D£\0\g), while the other continua are couple neither to 
the resonance nor radiatively to the ground state, and 
can therefore be entirely ignored. As a consequence, the 
transition amplitude ApE,g to a single final continuum /3 
through a multichannel intermediate continuum can be 
reduced to the coherent sum of two amplitudes: one for 
a single resonant intermediate channel, ApE,R,g, and one 
for a single non-resonant intermediate channel, Ai 3 E,D.g, 


Then, the amplitude becomes proportional to 

Af^ g oc w{za) ££ 2 exp (-r/2Ta), (84) 

where we made the approximation erfc(—iz) ^ 2 and 
we neglected the small term affST^. As expected, the 
transition amplitude decays exponentially with the time 
delay with half the lifetime of the resonant state. 

Notice that for negative time delays the resonant sig¬ 
nal still decays as the overlap of the pump and probe 
pulses (provided that the probe pulse is not itself in reso¬ 
nance with a transition from the ground state to a bound 
or autoionizing intermediate state). This latter circum¬ 
stance illustrates how, in a time-resolved formulation, the 
time ordering of photon exchange in the transition ma¬ 
trix elements translates to an actual order in the two- 
photon transitions, when the two photons belong to non¬ 
overlapping pulses. 

D. Multichannel case 

The results obtained thus far are valid only for the 
case of a single intermediate and a single final contin¬ 
uum channel. As long as the intermediate and final 
resonances are isolated, and if all the coupling matrix 


ApE,g = ApE,R,g + Ai3E,D,g- ( 88 ) 

A similar reasoning applies to the final states, since even 
in that case it is possible to identify a single final reso¬ 
nant continuum. However, due to the presence of multi¬ 
ple intermediate states (two different continua, the bound 
states and the autoionizing state), more than two decou¬ 
pled hnal continua can eventually be populated by means 
of a dipole transition. In any case, the final continua can 
be treated separately. If the final channel is not resolved, 
the individual contributions of all the final states to the 
total signal 'PE,g must be added incoherently, 

^Eg = '^\ApE,g\'^- (89) 

/3 

In other terms, the problem of several final channels can 
be treated as several problems of a single final channel 
(being it resonant or not). In conclusion, if all the rele¬ 
vant couplings with the intermediate and final resonant 
and decoupled unperturbed channels are available, the 
multichannel problem can be treated as a combination of 
the amplitudes given earlier in this section. 

Let us examine the case of one intermediate resonance 
in a multichannel continuum and no final resonances. Ac- 
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cording to the above and to Eq. (79), we can write 


1 


eiJa + 


- £Ea + * 


Cfia 


7 ) iQag-i)wiz^J 


(90) 


+ T^\t) A‘^-^Y.^PE.nOr,gW{z%\), 

n 

= T^\T)A-^^OpDODgW{zl^). (91) 

When taking the sum of the resonant and decoupled am¬ 
plitudes, the latter can be integrated in the first term in 
parenthesis of the former, giving rise to an effective com¬ 
plex q parameter. The overall amplitude, however, can¬ 
not be assimilated to a single resonant transition ampli¬ 
tude by simply redefining the parameters involved. Thus, 
in principle, the presence of a decoupled channel qual¬ 
itatively alters the finite-pulse resonant transition am¬ 
plitude. In the long-pulse limit, however, the situation 
changes, as the total transition amplitude becomes pro¬ 
portional to [compare with Eq. (82)] 




■pE,g OC roR + 


^ I „(21) 

e^ia + <7eff 


roR = 


OpdOdq 


(92) 


Efiia + « OpRORg 

The constant term trr expresses the strength of the 
dipolar coupling to the final continuum through the 
decoupled intermediate continuum \De) relative to the 
one through the (unperturbed) resonant continuum |i?e). 
Equation (92) can be rewritten as 


A 


21 

PE,g 


OC 


I (21) 

£gia + ^eff' 

^Eia + * 


(21) 1 + iroR 


( 21 ) 
9eff ■ 


(93) 


Therefore, in the long-pulse limit, the effect of multiple 
channels manifests itself as a simple modification of the 
effective complex q parameter, exactly as it happened in 
the case of a finite dipolar coupling between the interme¬ 
diate metastable state |a) and the final continuum j/JE), 
7o2 7^ 0. While in the latter case the modification of 
the effective q was different for the absorption and for 
the emission of the second photon, however, in the mul¬ 
tichannel case the variation of q is identical for the two 
paths. In principle, therefore, it is still possible to disen¬ 
tangle the two effects by comparing these two transition 
amplitudes. 


E. Multiple intermediate and final resonances 


The total transition amplitude (77) can be gener¬ 
alised to the case of several intermediate and final iso¬ 
lated resonances by adding to the common background 
term the individual contribution from the intermediate 
and final states plus the residual contributions from all 
intermediate-final resonance pairs. 
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iEa + * 
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eEb — i O0aOaE,g 
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21,ba Qag ^ 


PE,g 
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j{zf)-wizfJ 


€-Eb — * I ^Ea + 1 

+ ) [2* /3a + “ * + ^haisiab “ 0 “ Cba] \ ■ 


This approach has been employed to compute the spec¬ 
trum of a sideband comprising the 2p^ autoionizing 
state in the RABITT ionization of the helium atom from 
the ground state, when both the lower and the upper 
harmonics contributing to the resonant sideband were 
themselves in resonance with the SP 22 the 
states, respectively [10]. 


F. Multiphoton transitions 


The n—th order finite-pulse transition amplitude (7) is 


.{n) ^ ^ 

3® (27r)v“i 


Sic,fg-n'„)l[[Fiujl)M 

2=1 


X 


n—1 

x{f\Ol[[G+{Eg+n'^0]\g), (94) 

i=l 


where fl' = X]j=i 4 ^.nd the factors in the last oper¬ 
ator product are assumed to be ordered from right to 
left. As long as the on-shell approximation is justified, 
the techniques employed in Sec. Ill B for the two-photon 
transition matrix element can be used also to compute 
n—th order transition matrix element. Furthermore, if 
no more than one intermediate resonance contributes to 
the transition, the procedure followed in Sec. IIIC can 
be subsequently applied to evaluate the folding with the 
field. A particularly relevant example that meets these 
conditions is the absorption of one pump photon uii fol¬ 
lowed by that of two probe photons i 02 , when only the 
first intermediate continuum \tjjae) is resonant while the 
second intermediate continuum |/3e) and the last contin¬ 
uum I ye) are not. In this case, the three-photon transi¬ 
tion matrix element is 

= mOGUEg+uj[+uj'2)OGUEg+uj[)0\g) 

_ _ ^ _ yw(21) (, / ) 

- Eg+uj[+oj'2-E + t0+^f^Eg^'^^)- 


(95) 
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The transition amplitude then becomes 


^221 _ ^^ 7/3 




duj. 


iO^(^ 
W2 . 


, F2iuJEg - - 0^2)-^2 K) 

Eg + Uj[ + UJ2 — E + iO + 
(( 21 ) 


dujM)^E g(^)Ei{uj)E^{uJEg - uj),{9&) 


where in the last passage we have assumed that the spec¬ 
trum of the absorption component of the probe pulse is 
localised around 102 and we used the convolution theo¬ 
rem. This means that the three-photon amplitude g 
is equal, apart for a multiplicative factor, to the two- 
photon amplitude in which the frequency and spectral 
width of the probe field are larger by a factor of 2 and 
y/2, respectively. 


IV. RESONANT RABITT SPECTRUM OE 
HELIUM 

In this section we illustrate the capabilities of the 
finite-pulse two-photon resonant model by computing the 
RABITT photoionization spectrum of the helium atom 
from the ground state to the energy region be¬ 

tween 30 eV and 40 eV above the first ionization thresh¬ 
old, which features the series of metastable doubly ex¬ 
cited states that converge to the N = 2 threshold. Fig¬ 
ure 9 shows the energy levels of helium in the region of 
interest and illustrates schematically the radiative cou¬ 
plings that must be plugged in the model to reproduce 
the RABITT spectrum of the atom when the harmonics 
can be resonant with the first two ^P° bright autoion- 
izing states, and sidebands can populate the final 2p^ 
state. Helium is an ideal candidate to investigate 
atomic transitions through autoionizing states because it 
is amenable to an accurate ab initio description in the 
presence of external light pulses. Furthermore, the ion¬ 
ization continuum of helium has been the subject of in¬ 
tense study for more than fifty years [56]. In particular, 
the N = 2^V° and series of autoionizing states have 
been investigated both experimentally [51, 57-60] and 
theoretically [61-71], and many of their properties, such 
as positions, width, and q parameter from the ground 
state, are well known. 

In the intermediate states of the model we included 
one ^P° intermediate continuum, IsEp, with the two iso¬ 
lated resonances spj and sp^ [62] , and one intermediate 
bound state, ls2p, which contributes significantly to the 
excitation amplitude of the final 2p^ state, owing to 
the strong dipolar coupling between the Is and 2p or¬ 
bitals. In the final states of the model we included the 
IsEg continuum, featuring the 2p^ metastable state, 
and the IsEd continuum. For the latter, we did not 
include any resonance, as the most relevant one, also with 
dominant configuration 2p^, lies very close in energy to 
the SP 2 state and hence it is not reached by any sideband 


IsSa IsEp IsEd 



FIG. 9. (Color online). Scheme of the essential states involved 
in the RABITT ionization of the helium atom in the region 
of the N = 2 autoionizing states, together with the relevant 
radiative couplings between them that must be taken into ac¬ 
count to reproduce the pump-probe photoelectron spectrum 
with the finite-pulse resonant two-photon model described in 
the text. The dashed lines indicate the non-radiative coupling 
between the resonant and the degenerate continuum states. 


within the chosen range of IR frequencies. The sp~ and 
2pnd states, as well as higher terms in the sp+ series, 
which are all narrow and have a small dipole coupling 
with the ground state if compared with the spy^ states, 
are not expected to affect significantly the sideband spec¬ 
trum near or below the 2p^ state and were therefore 
not included in the model. 

The position, width and q parameter from the ground 
state of the two intermediate resonances, as well as the 
background photoionization cross section, can be taken 
from the literature, where one can find also the posi¬ 
tion and width of the final state and the energy of 
the bound ls2p state. Even with these data, there are 
still 14 independent parameters not reported in the lit¬ 
erature that are in principle needed to apply the model: 
the two continuum-continuum couplings, OisEi, 1 sEj,] the 
relative strength of the direct dipole coupling of the two 
intermediate autoionizing states with the two final con- 
tinua, PisEt,sp+ (4 parameters); the dipolar 

coupling to the two final continua through the interme¬ 
diate ls2p bound state, OisE^,is2pC’is2p,is2; the q param¬ 
eter for the excitation of the final resonance from the 
intermediate bound and the two intermediate autoion¬ 
izing states, 92p".u2p’ ^ 2 p\sp+^ ^ 2 p\sp+^ and, conversely, 
for the excitation of the spy^ intermediate resonances 
from the 2p‘^ final metastable state, 9^~p+2p2, Qs~p+, 2 p^'^ 
the dipole coupling between the final metastable state 
and the intermediate continuum, C*2p2 The value of 
the residual parameters ^2p=.sp+> ^ 2 p^,sp+^ ^ 2 p^,sp+y and 
(^ 2 p 2 sp+, can be determined from the previous ones with 
















17 


the additional assumption O„o — 0„~ 2 + , which 

'‘P i'5P2/3 ^P ’®P2/3 

is justified by the strong dipole coupling between dou¬ 
bly excited states compared to that between an = 2 
doubly excited state and a IsE^ continuum (the latter 
optical transition, being itself a double excitation, is pro¬ 
hibited within the quasi-particle approximation). In the 

same spirit, we can assume 0-+ ,0 = 0 + with 

*P2/3’ P *P2/3’ P 

which q - + , 2 become derived quantities, thus reduc- 

ing to 12 the total number of independent parameters. 
Finally, we assume ^ 1, so that the product 

Oisee,is 2 pOis 2 p,is 2 q 2 p 2 comes as a single parameter 
"l 2 V'(-is 2 thus reducing the total number of param¬ 

eters to fl. 

In one of our previous works [48] , we developed a soft- 
photon model for the non-resonant ionization of helium 
in which the continuum singlet states of the atom were 
approximated by the product of the Is ground state of 
the He^ parent-ion with free spherical waves, 

aP - Pal+'Pi2 ■ n ^ ^ 

(pis, Elm = —- [ri)y —je{kr2)Yem[r2), 

(97) 

where jg are spherical Bessel functions, Yim are spherical 
harmonics, and E = fc^/2, while the ground state was 
approximated by the Is^ configuration, where, following 
Slater’s prescription, the Is orbital was an hydrogenic 
wavefunction with effective charge Z = 1.7, 

p /ov5/2 

With this model we were able to predict, with quantita¬ 
tive accuracy, the background distribution of the photo¬ 
electrons generated by the interaction of the atom with 
sequences of XUV-pump and IR-probe pulses, even in 
the presence of autoionizing resonances and for large 
IR intensities, provided that the whole spectrum (which 
occupied a limited energy region approximately 1 a.u. 
above the Is threshold) was scaled by a constant factor 
C, of the order of unity, which accounts for the known 
difference between the hydrogenic ionization cross sec¬ 
tion and the one predicted by the first Born approx¬ 
imation. Since we are currently considering the same 
region of the photoelectron spectrum, it is justified to 
estimate the continuum-continuum couplings OisEe,isEj, 
within the soft-photon approximation, further reducing 
the total number of free parameters to 9. 

To determine the values of these parameters, and to 
subsequently verify the prediction of the model against 
reliable benchmarks, we carried out accurate ab ini¬ 
tio simulations based on the numerical solution of the 
time-dependent Schrbdinger equation [10, 32-34] where 
the wave function for the two-active-electron system 
is represented in a B-spline [72, 73] bipolar-spherical- 
harmonics [74] close-coupling basis [75-78], for both the 
time propagation and the asymptotic analysis of the 
single-ionization wave packet [32, 79, 80], while the time- 
step propagation is carried out with a Krylov unitary 



photoelectron energy (eV) 


FIG. 10. (Color online) Pliotoelectron spectra obtained ab 
initio (thick gray solid line in the background) and with the 
model (thin black line in the foreground) for four different 
pulse sequences useful to calibrate the model parameters, a) 
non-resonant and b) resonant one-photon photoelectron spec¬ 
trum obtained with using the XUV APT alone. The first spec¬ 
trum (a) was used to determine the global scaling factor C, 
while spectra like the one in (b) were used to confirm the pa¬ 
rameters of the two ^P° autoionizing state, c) from the lower 
sideband of the spj resonance generated with disjoint pump 
and probe pulses, we could estimate the resonance-continuum 
P radiative couplings, d) from the resonant shape of the sp^ 
state upper sideband, which strikes the 2p^ state, we could 
determine the resonance-resonance dipolar coupling. 


approximation [81] to a second-order split exponential 
time-evolution operator in velocity gauge, and imple¬ 
mented in a parallel code which makes use of the PETSc 
numerical library [82-84]. In the simulations, both the 
XUV APT and the IR probe have a duration (fwhm) of 
6 fs, while their peak intensity is Ixuv = 5 GW/cm^ and 
/m = 10 GW/cm^, respectively. The individual attosec- 
ond pulses in the train have central energy Wxuv = 57 eV 
and a duration of 250 as; consecutive pulses are separated 
by half the nominal IR period, Tm = 27r/a;iR. As we are 
particularly interested in the transition that, through the 
sp 2 and sp^ (^P°) DESs, populates the optically for¬ 
bidden 2p^ (^S'®) DES, we performed simulations for IR 
frequencies ranging from Wia = 1.455 eV to Wia = 1.485 
eV. 

Using the XUV APT alone, we could determine the 
value of the scaling constant C to match the non-resonant 
background of the model to that of the ab initio pre¬ 
diction (which agrees with the absolute value of the 
background photoionization cross section reported in the 
literature). With C = 1.2, the model and ab initio 
backgrounds are in excellent agreement across an en¬ 
ergy domain of ^ 4wia (see Pig. 10a). Prom the one- 
photon spectrum in the energy region where the 
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and resonances are present, we determined position, 
width and q parameter for the two autoionizing states: 
= -18.86 eV, r + = 0.037 eV, g ~+ = -2.77, 

4 + = -15.34 eV, = 0.0082 eV, <? ~+ = -2.58 

(see Fig. 10b), in agreement with the values reported in 
the literature [71]. 

We obtain the parameters d, ^ + from the ab ini- 

tio background spectrum, resolved with respect to the 
orbital angular momentum £, of the two intermediate 
resonances for a time delay at which the pump and the 
probe pulses do not overlap since, as we already saw in 
the Sec. IllCOj, in this way the homogeneous contribu¬ 
tion of the intermediate state vanishes. Notice that one 
could retrieve the same parameters also from the exper¬ 
iment by measuring the photoelectron spectrum at two 
different ejection angles. We determine both f3isEp,2p^ 
and 72 p^<_is 2 p<-is 2 matching the parameters of the 
asymmetric resonant profile in the two-photon excitation 
of 2p^ with a pair of overlapping pump and probe pulses 
in which the harmonics are tuned out of resonance with 
respect to the intermediate autoionizing states. Finally, 
to determine q + „ 2 , we look at the 2 p^ resonant pro- 

^P2/3^‘^9 

file in the sideband of the ab initio spectrum for non¬ 
overlapping APT and IR probe pulses, where alterna¬ 
tively the lower and the upper harmonics are in resonance 
with the SP 2 and the sp^ state, respectively. 

In Table 1 we report the full list of the parameters that 
gave the best match with the ab initio spectra discussed 
above. 

TABLE I. Radiative parameters for the two-photon resonant tran¬ 
sitions model in helium. Atomics units are used. 



Ea 

r. 

<?a 

PlsEg,a 

PlaE^,a 

a=sp^ 

-0.693 

1.37[-3] 

-2.77 

-0.003 

-0.003 

a=spj 

-0.564 

3.01 [-4] 

-2.58 

-0.003 

-0.01 



Et, 

Ft 

%,aj>+ %,sp+ 0b,l3Ep 76<ls2p<-ls2 

b=2p2 

-0.622 

2.16 [-4] 

153 20 -0.003 4255 


Notice that almost all the parameters of the model 
were determined independently of each other by com¬ 
parison with a minimal number of well defined selected 
numerical experiments. Alternatively, the parameters 
could have also been determined by comparing with ac¬ 
tual time-unresolved experiments [87]. Once these val¬ 
ues are determined, the model can reproduce the photo¬ 
electron spectrum for several values of the IR frequency 
and pump-probe delay in the general case of partly over¬ 
lapping pulses with multiple intermediate and one final 
resonant states, with no residual freedom to adjust the 
outcome. 


The left panel of Fig. II shows both the ab initio (left 
panel) and the model (central panel) prediction of the 
photoelectron spectrum as a function of the pump probe 
time delay for the sidebands 8638-42 at a fixed IR fre¬ 
quency of = 1.466eV. Positive time delays indicate 
the XUV pulse train comes first. At the mentioned fre¬ 
quency, the 4Ist harmonic is resonant with the spj i}P°) 
doubly excited state {E^p+ « 35.55eV). The intermediate 
resonance induces a local phaseshift, in opposite direc¬ 
tions, of the resonantly populated sidebands SB 40 and 
SB 42 , an effect that is not present when a sideband is 
populated via non-resonant paths only, as in the case 
of SB 38 . Model and ab initio spectra look essentially the 
same. To illustrate qualitatively how the model compares 
with the ab initio theory, on the right panel of Fig. II we 
show the two predictions for the energy integrated spec¬ 
trum of the two sidebands as a function of the time delay, 
and they are indeed found to be in excellent agreement. 
The vertical white lines in the two first panels denote the 
maximum of the sideband signal and show how the ab¬ 
solute value of the local phaseshift 5ip{loi-r,t) of the two 
resonant sidebands increases with the time delay, 

Isb{t) oc cos{2wiar -f (5v?(Wm,T)} . (99) 

This means that the resonance introduces a modulation 
of the RABITT beating frequency itself. In this scenario, 
the concept of a global RABITT phase loses its meaning. 
In the cases we examined, however, the phase deviation is 
well approximated by a linear interpolation, (I<p(a;iR,T) Ri 
(5(/3o(wm) -I- so 

IsB OC COS {[2wia -|- (5a;(wia)] t -|- (5<po} ■ (100) 

The local phaseshift is affected by the apparent phase- 
shift at r = 0 as well as by the modulation of the fre¬ 
quency, ^aj(a;ia). As discussed in Sec. II Oe, the mod¬ 
ulation of RABITT beating frequency appears even in 
absence of intermediate resonances, as a result of using 
finite pulses. The latter non-resonant effect, however, is 
always a shift towards the red, it does not depend much 
on the IR carrier frequency, and it rapidly disappears 
as longer pulses are employed. The resonant modula¬ 
tion of the sideband frequency, on the other hand, in¬ 
duces opposite shifts in the two resonant sidebands, it 
depends strongly on the detuning of the resonant har¬ 
monics from the intermediate autoionizing states, and 
it becomes sharper when longer pulses are used. The 
model prediction for the frequency modulation, obtained 
by Fourier analyzing the energy integrated sideband sig¬ 
nal, is 5uj = —0.073 eV. The non-resonant redshift asso¬ 
ciated to the use of a 800 nm 6 fs probe pulse is com¬ 
paratively large, Sujnr = —0.038 eV. Both the total and 
the non-resonant values are in agreement with those from 
the ab initio calculation (the latter being estimated from 
the non-resonant sideband SBss). By taking the differ¬ 
ence between the total and the non-resonant values, the 
resonant contribution to the sideband frequency modu¬ 
lation, due to the spj doubly excited state, is estimated 
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FIG. 11. (Color online). Left panel: photoelectron spectrum as a function of the pump-probe time delay for sidebands SBss, 
SB 40 and SB 42 of the driving frequency cjir = 1.467 eV. Right panel: energy-integrated photoelectron spectrum as a function 
of the pump-probe time delay for both model and ab initio simulations. The presence of the sp^ DES (~ 35.5 eV), which is 
populated by the 41st harmonic (not shown) shifts sidebands SB 40 and SB 42 in opposite directions. Furthermore, due to the 
finiteness of the pulses used, the resonance induces a frequency modulation which can be seen by comparing the separation 
between the white lines that indicate the maxima of the sideband oscillations. 


as Sui^p+ = —0.035 eV, which corresponds to a change in 
the RABBIT period of 17 as. 

Figure (12) shows the photoelectron spectrum of side¬ 
bands 40 and 42, as a function of the IR carrier frequency, 
for five different pump-probe time delays. The agreement 
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FIG. 12. (Color online). Photoelectron spectrum as a func¬ 
tion of the driving laser frequency for SB 40 and SB 42 at five 
different stages of the sideband oscillation (r = Tr corre¬ 
sponds to the minimum). Upper panels show the ab initio 
calculations and lower panels show the model results. 



between ab initio and model, again, is excellent. In par¬ 
ticular, the sideband resonantly populated from below 
(SB 42 ) shows a maximum to the left and a minimum to 
the right of the central resonance frequency {huj Ri 1.466 


eV), while the opposite is true for the sideband that is 
populated from above. This feature is responsible for 
the apparent phaseshifts for the two sidebands, which are 
shown in the first two panels of Fig. 13 and were obtained 
by Fourier analyzing the spectrum in the time-delay in¬ 
terval T G [0 ,T/r/ 2]. The sp 2 and spg, populated by 
H41 and H43, respectively, give rise to resonant struc¬ 
tures in the apparent phaseshift that are located at IR 
frequency close to the resonance condition of the each 
DES with the corresponding harmonics. As discussed 
in Sec. Ill, the overall phase excursion depends on the 
parameters of both the resonance and the pulses used. 
In the present case, the larger dipole matrix element of 
SP 2 with the ground state, compared with that of sp'^, 
makes the former dominate the shape of the profile, al¬ 
though the peak for the sp^ DES can also be recognised. 
Einally, in the last panel of Eig. 13 we show the phase 
of the beating of the H 39 integrated harmonic signal. In 
this case, the resonance profile arises from the interfer¬ 
ence between the direct one-photon ionization amplitude 
from the ground state and the three-photon amplitude 
for the resonant absorption, from the ground state, of 
one XUV photon of the H 41 harmonic followed by the 
stimulated non-resonant emission of two IR photons. As 
we discussed at the end of Sec. Ill, this latter amplitude 
can be easily computed with a straightforward multipho¬ 
ton extension of the finite-pulse resonant model. The 
good agreement with the full-fledged ah initio simulation 
for this process certifies that such extended model works. 
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FIG. 13. (Color online) Apparent phaseshift of the integrated- 
signal beating of sidebands SB 42 (left panel), SB 40 (central 
panel), and harmonic H 41 (right panel), with respect to the 
non-resonant SBas sideband phase, as a function of the fun¬ 
damental IR carrier frequency. In all three cases, the largest 
phase excursion occurs when the intermediate harmonic H 41 is 
resonant with the spj DES. In SB 42 , the effect of through 
H 43 is visible as well. The resonant profile of H 41 , which is 
comparable to the one of SB 40 , results from the interference 
between the direct one-photon excitation amplitude of the 
continuum from the ground state with the three-photon am¬ 
plitude, in which the resonant absorption of one pump pho¬ 
ton from H 41 is followed by the stimulated emission of two IR 
probe photons. 


nant transitions in any atoms, molecules or solids sus¬ 
ceptible of a description in terms of a finite number of 
free-particle channels and metastable states. Think, for 
example, of the radiative excitation of image-potential 
states on metal surfaces, which, on the one side, can de¬ 
cay by tunnelling to the conduction-band, and on the 
other side, can exchange a further photon and be liber¬ 
ated to either the metal or to the vacuum (photoemission 
channel) [85, 86]. 

When used as a phenomenological tool, the model can 
be employed to extrapolate, from time-resolved exper¬ 
iments with table-top laser apparatuses, the radiative- 
coupling strength between short-lived excited states, such 
as autoionizing states in heavier rare gases, which can 
be hard to obtain otherwise, either theoretically, due 
to the challenging role of electronic correlation, or ex¬ 
perimentally, due to the need of coupling lasers to a 
synchrotron x-ray beamline [87]. Conversely, once the 
model is parametrised, it can be used as a computation¬ 
ally inexpensive alternative to the numerical integration 
of the TDSE. This is for example the case of photoe¬ 
mission studies conducted with the long, coherent and 
intense XUV pulses that became recently available at 
seeded XFEL facilities [88-90]. 



V. CONCLUSIONS 


Appendix A: Faddeeva function. 


In this work we have presented a new finite-pulse 
model for resonant two-photon transitions, which can 
be applied to simulate, at a negligible computational 
cost, attosecond pump-probe photoionization processes 
in atoms. The model, which extends the one presented 
in [10], accounts for multiple intermediate and final chan¬ 
nels, as well as the possible presence of multiple isolated 
resonances. Further generalisation to higher order tran¬ 
sitions has been outlined. 

We used the model to explain the physical origin of res¬ 
onant phase profiles in two-photon ionization spectra as a 
function of the pump-resonance detuning. In particular, 
we showed that, if the intermediate states comprise a sin¬ 
gle continuum and a resonance not radiatively coupled to 
the final continuum states, the phase of the two-photon 
amplitude coincides with that of Fano one-photon transi¬ 
tion, while in presence of multiple intermediate continua 
or of a finite radiative coupling between the intermedi¬ 
ate resonance and the final continuum states, the phase 
experiences a continuous excursion with a net variation 
that can be either 0 or 27r. Furthermore, we showed that, 
when very short pulses are used, the beating frequency of 
the sidebands in the RABITT pump-probe scheme dif¬ 
fers from the nominal 2ujyb_ value. The results obtained 
with the model are found to be in quantitative agree¬ 
ment with virtually exact ab initio simulations for the 
RABITT photoionization of helium in the region of the 
N = 2 doubly excited states. 

Even if benchmarked here against helium, the model 
is applicable to study time-resolved multiphoton reso- 


In this appendix we derive the general analytical ex¬ 
pression for the two-photon transition amplitude between 
an initial state \g), with energy Eg, and a final state \/3E), 
with energy E, due to the absorption/emission of a pho¬ 
ton from a first Gaussian pulse Ei, centered in t\ = 0, 
followed by the absorption/emission of a photon from a 
second Gaussian pulse E 2 , centered in t 2 = ti + t = t, 


•^^PE,g — 


da;F2(wEg - a;;r)Fi(w)AI/3E,g(a;), (Al) 


under the hypothesis, recurrent in the derivation of the 
model illustrated in Sec. Ill, that the two-photon matrix 
element AI/ 3 _E,g(w) has an isolated simple pole at w = 
Ea — Eg, where Ea £ C, ImFo < 0, and that, to a very 
good approximation, (bj + Eg — Ea)AipE,giy>) is constant 
in the region where the product E 2 {uJEg — i^', 'r)F'i(w) does 
not vanish. 


MpE,g{u}) 


^PE,g 

UJ -UJag' 


(A2) 


In the following, to consider all the possible cases at 
once, we will indicate both the absorption and the 
emission spectral components 29 of the n-th Gaussian 
wavepacket (28) with the single expression 




An 


2(7 ni 


g-iV3sgn(w„)gii.;t„g' 



(A3) 


where the absorption/emission components are differen¬ 
tiated by attributing to uq a positive or negative sign. 
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respectively. The two-photon transition amplitude (Al), 
thus, becomes 


J21 _ 

■^l3E,g - 


— sgn(a; 2 ) p-ivi sgn(wi) 


-*7:— e 
2(72 


2 (Ti 


(A4) 


X T, 


PE,g 


dw-e ^^2 e 


^ag 


After some lengthy but straightforward algebraic pas¬ 
sages, it is possible to cast the previous result in the 
following form 


x (A5) 

^ 2(J2 2(Ji 


X exp 


2(7^ 2(7, 


T" .(72 T (5 

- 2 ~ *-f ) X 


(7l (7t (7 


/ OO 

doj 

-OO 


exp 




Wi -I- W - Wag 


where we have introduced a convoluted spectral width 
(7 = a/( 7 ^ -I- ( 7 | and temporal width at = \J erf ^ -I- a^^ 
(notice that a = ( 7 i( 72 ( 7 t), as well as the nominal detun¬ 
ing 5 = Eg + LOi + UJ 2 — E- By performing the change 
of variable x = + fjf + *^), the integral in 

Eq. (A5) can be expressed in terms of the Faddeeva spe¬ 
cial function w{z) = e~^ erfc(—iz), which, in the upper 
half of the complex plane, admits the following integral 
representation (see §7.1.3-4 in [55]), 


w{z) = — - dt, lm[z] > 0. (A 6 ) 

E J—OO ^ t 


Indeed, by introducing the dimensionless complex vari¬ 
able 


21 ^ 

“ V2 


\U-is- 

■a) 

^ai 

LV 0-2 




(A7) 


the integral in Eq. (A5) can be expressed as 


.OO exp 
/ duj - 

*/ —OO 




- (7tW -I- 


(7i S 


iJl CJ — LJa 


= -mw{za). 


(AS) 

Notice that to establish the correspondence between the 
integral in Eq. (A5) and the r.h.s. of (A 6 ), one must 
continuously deform the integration path from the initial 
real axis to a final re-defined real axis without crossing 
the pole, which requires t < at ImEa- Once the integral 
is written in terms of the Faddeeva function, however, the 
expression is valid for any value of the time delay, since 
the Faddeeva function is defined on the whole complex 
plane by analytic continuation. The time-ordered two- 
photon transition amplitude finally becomes 


a21 _ A 1 A 2 -iya sgn(a; 2 )p-iyi sgn(a;i) ^ iuJ 2 ^ 

■^ISE,g- ® ^PE,ge 


X exp 


-b 
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